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Abstract 

We revisit the globally coupled map lattice (GCML). We show that in the 
so called turbulent regime various periodic cluster attractor states are formed 
even though the coupling between the maps are very small relative to the non- 
linearity in the element maps. Most outstanding is a maximally symmetric 
three cluster attractor in period three motion (MSCA) due to the foliation of 
the period three window of the element logistic maps. An analytic approach 
is proposed which explains successfully the systematics of various periodic- 
ity manifestations in the turbulent regime. The linear stability of the period 
three cluster attractors is investigated. 
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I. INTRODUCTION 



Recently there has been much progress in the study of synchronization of nonhnear maps 
and flows [pHT^. This may lead to the clarification of the intelligence activity supposed 
to come from the synchronization among the neurons in the neural network. Especially the 
globally coupled map lattice (GCML) may be considered as one of the basic models for the 
network systems expressing their characteristic limits. In its simplest form, all elements 
interact among themselves via their mean field all to all with a common coupling, and each 
of the element is a simple logistic map with a given nonlinearity. It may be regarded as a 
natural extension of the spin-glass theories |T^JT^ to the nonlinear dynamics. Even though 
the simplest GCML has only two model parameters — the common nonlinearity parameter a 
and the overall coupling e — it exhibits a rich variety of interesting phases on the parameter 
space corresponding to various forms of synchronization among the maps determined by the 
balance between the randomness specified by a and the coherence by e. 

In this article we revisit the turbulent regime of GCML, which is a regime in the pa- 
rameter space with high a and very small e. The main interest in this regime has been so 
far focused to the so called hidden coherence 0,^ . It is a phenomenon that the fluctuation 
of the mean field of the maps in evolution does not cease at large system size. The mean 
field distribution obeys the central limit theorem (CLT) but not the law of large number 
(LLN) We show that the dynamics of GCML in this regime is a foliation of that of the 
element logistic maps and that various periodic cluster attractors are formed even though 
the coupling between the maps is set very small. We show that the regions which may be 
described by the hidden coherence do exist but that it is a very limited part of the parameter 
space. 

We organize our discussion in three parts. First, we present the results of an extensive 
phenomenological survey of this regime and list evidences of periodicity manifestations due 
to the periodic windows of the element logistic map. Most outstanding is the onset of 
period three cluster attractors. The turbulent regime is, if we may say, a bizarre region with 
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many faces — drastic periodicity manifestations as well as almost perfect randomness under 
the hidden coherence. At the periodic or quasi-periodic attractors, the mean field evolves 
controlled by the scale of the cluster orbits and the LLN is naturally violated. We also 
present a remarkable data which shows that the GCML at the large system size acquires a 
high sensitivity to the periodic windows of the element map. Second, we present an analytical 
approach which successfully explains the systematics of the periodicity manifestations. We 
present a tuning condition which limits the system parameters with which GCML cluster 
attractor states of a given periodicity may be formed. The key to obtain the condition is 
the introduction of the maximally symmetric cluster attractor (MSCA), which is a solution 
of minimum fluctuation in the mean field. It corresponds to the known state of two clusters 
in opposite phase oscillation which is formed in the ordered phase of GCML We verify 
the validity of the condition in detail and show that the foliation is the governing dynamics 
of this regime. 

Third, we show that the period three cluster attractors formed in the turbulent regime are 
linear stable and investigate how their stability changes by the coupling e and the population 
ratios. We in particular derive algebraically the e value for the formation of the most linearly 
stable bifurcated-MSCA (MSCA*). 

The organization of this article is as follows. In Sec. II we briefly review various GCML 
phases and locate the turbulent regime in the parameter space. We then summarize the 
known facts of this regime. No originality is claimed here. We then briefly compare them 
with our results. In Section III we present our phenomenological findings including the 
period three MSCA and associated fewer cluster attractor. In Sec. IV we present an analytic 
approach which explains successfully how the periodic windows of the element map control 
the GCML dynamics in the turbulent regime. In Sec. V we investigate the stability of the 
cluster attractors. We conclude in Sec. VI. 
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II. GCML PHASE STRUCTURE AND THE TURBULENT REGIME: A REVIEW 

We study in this article the simplest GCML which is a system of maps evolving by 
x,{n + l) = {l-e)f{xi{n))+€h{n), = 1, ■ ■ ■ , iV), (2.1) 
and the mean field h{n) of maps is defined as 

-I N 1 TV 

= T7E/(^^W) = ^T.^^ir^ + l). (2.2) 
1=1 1=1 

In the first step, all Xj are simultaneously mapped by a nonlinear function /. The function 
/ could be distinct maps (heterogeneous GCML) but in this article we consider the simplest 
case that / is a logistic map f{x) = 1 — ax"^ common to all variables (homogeneous GCML). 
The nonlinearity of / generally magnifies the variance among the maps. The larger the 
parameter a is, the more strongly the variance is enhanced. In the second, the maps undergo 
interaction between themselves with a global coupling constant e. Here every /(xj) is pulled 
to their meanfield h{n) at a fixed rate 1 — e. The larger the coupling e is, the more strongly 
the maps are driven into synchronization. 

The model is endowed with various interesting phases under a subtle balance between 
the two confiicting tendencies. The phase diagram in the a, e plane was explored by Kaneko 
|T|]. Let us explain the phases choosing a = 1.80 for definiteness. This is far above the 
criticality a = 1.401 ■ ■ ■ to the chaos for a single logistic map. (i) For a sufficiently large e 
(> 0.38) the maps are strongly bunched together in a cluster in the final attractor and evolve 
chaotically as a single logistic map. This is the coherent phase, (ii) For e = 0.22 — 0.30 
the interaction via the mean field can no longer exerts strong bunching and the final maps 
divide into two clusters. The maps in each cluster are still tightly synchronizing each other, 
and the two clusters mutually oscillate opposite in phase. This phase turns out as a solution 
of a minimum fiuctuation in the mean field and called as the two-clustered ordered phase. 
(iii) For smaller e, the number of final clusters increases but it remains independent from 
the total number of maps in the system. The typical number of clusters at various e ranges 
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is indicated in the phase diagram |T|]. (iv) Finally, for very small e the number of clusters 
is in general proportional to A^. This region is called as the turbulent regime. This is the 
target region of our analysis. 

It is known that in the turbulent regime maps evolve almost randomly at small lattice 
size N and that there occurs a subtle correlation — a hidden coherence — at large N . But 
as we show below there actually emerge drastic global periodic motions of maps if the e 
takes certain values for a given a. Let us first briefly review previous observations in the 
literature. 

(i) The final state of GCML in this regime iterated from a random configuration consists 
of maps and tiny clusters, each moving chaotically due to high nonlinearity. The number of 
elements (maps and clusters) is proportional to the number of whole maps in sharp contrast 
to the ordered regime 0. 

(ii) There emerges certain coherence between elements when the size is large P,^. If 
Xi{n) {i = 1, ■ ■ ■ A^) are really independent random variables following a common probability 
distribution, the mean squared deviation (MSD) of the mean field h{n) {5h^ = {h"^) — {h)'^ 
with (■ ■ ■) here meaning the long time average) should decrease proportionally to by the 
law of large numbers (LLN) and the h{n) distribution must be a gaussian for sufficiently large 

by the central limit theorem (CLT). However there is a certain threshold in A^ (depending 
on both a and e) above which MSD ceases to decrease even though the distribution remains 
gaussian; CLT holds but not LLN 0. This reflects some hidden coherence between the 
maps in evolution. In fact LLN is restored when a tiny noise term is introduced in each map 
independently . 

(iii) The violation of LLN reflects that the map probability distribution p{x) depends 
on time. Indeed a noise intensity analysis of ensembles successfully proves LLN 0. If 
LLN should hold in the time average, p{x) would have to be a fixed point distribution of 
the Frobenius-Perron(FP) evolution equation ||T5|. It has been argued that the fixed point 
distribution may be unstable due to the periodic windows of the logistic map though 
this point is controversial. For instance, on tent maps, the same instability occurs but no 
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periodic windows are present |T6|. The coherence manifests itself in the mutual information 
0. On the other hand the temporal correlation function similar to the Edwards- Anderson 
order parameter for the spin glass ||14| decays to zero exponentially. Thus it may not be due 



to freezing between two elements 0. 

The hidden coherence was found in the statistical analysis of the mean field fluctuation 
0. But there has been no report of an extensive statistical analysis which covers the whole 
turbulent regime as well as a wide range of the system size. And once we have done it, we 
are faced with a bizarre feature of the turbulent regime; the hidden coherence is one thing 
but there occur also drastic global periodicity manifestations. The above lists are correct 
but need reservations. 

As for (i) there is a need for a careful reservation on the coupling values. We show below 
that, when the coupling e takes small but tuned values for a given a, the maps again - like in 
the ordered phase - may split into a few bound clusters in periodic motion. The most striking 
manifestation of periodicity in this form is the states of almost equally populated clusters 
mutually oscillating in the same period with the number of clusters. We call this type of a 
periodicity manifestation as a maximally symmetric cluster attractor (MSCA) and present 
below the tuning condition for it. We label such a cluster attractor by the periodicity and 
the number of the clusters. For instance, we call the outstanding period three symmetric 
cluster state as p3c3 MSCA. There also occurs bifurcated p3c3 MSCA. The MSD of the h{n) 
distribution is very small in the MSCA or its bifurcated state because of the good population 
symmetry among the clusters. At slightly larger e, we observe that the number of clusters 
decreases while the orbits are approximately kept. The cluster attractor of this type {p > c) 
leads to large MSD, which is independent from N. 

As for (ii) we show below not only LLN but even CLT is violated in almost all regions in 
the turbulent regime. We pin down the very limited regions in the turbulent regime where 
the CLT holds with violated LLN; only there the term hidden coherence may be used. 

As for (iii) the decay exponent of the temporal correlator of the mean field fluctuations 
gradually decreases with the deviation of the coupling from the tuned value. Accordingly 
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the h{n) distribution successively changes its shape from the highest rank sharp delta-peaks 
down to the MSD-enhanced Gaussian distribution — the hidden coherence. This indicates 
that the hidden coherence at the MSD valley may be the most modest periodicity remnant, 
being elusive due to high mixing. 

The GCML can be defined in a one line equation but its turbulent regime challenges 
us with so many faces ranging from a manifest periodicity to the hidden coherence. We 
consider that it is important to explore the systematics of periodicity manifestations by 
an extensive statistical survey and present a sorted list of phenomenological observations. 
Below we firstly devote ourselves to this task. 

III. PHENOMENOLOGY OF THE TURBULENT REGIME 

A. Systematics in the mean field fiuctuations 

We start with an analysis of the distribution of the mean field fiuctuation in time. In 
Fig. |l] we show its MSD at a = 1.90 as a surface over the N — e plane, which overlays a 
density plot of the periodicity-rank of the h{n) distribution. In order to set sufficiently fine 
grids for the surface, the system size is limited in the range < 4 x 10^. For a wider range 
of A^, we show in Fig. ^ the sections of the MSD surface at in powers of ten up to 10^. 

1. Peak-valley structure of the MSD surface 

First let us discuss the MSD surface and its sections. The linear edge of the surface at 
e = is of course due to LLN. For a non-zero but very small e 0.01), the LLN still 
holds to a good approximation but otherwise we can clearly see that the surface has many 
peaks along the A^— axis — the violation of LLN in the time-seriesQ. There is a prominent 

^This does not imply the real violation of the LLN in the ensemble average The violation of 
LLN means here that there is a larger fluctuation in the time series of h{n) than expected by LLN. 
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peak at £ ~ 0.040 — 0.050 — an extreme violation of LLN — and in this e range the MSD is 
in excess even for ^ 10^. In front of the peak there is a deep valley around e ~ 0.035. 
We show below these peak and valley are respectively induced by p3c2 cluster attractor and 
p3c3 MSCA (and its bifurcated state). Apart from these, the MSD surface systematically 
shows successive peak-valley structure at large A^. We can see clearly in Fig. ^ how this 
structure is formed with the increase of A^. At MSD peaks the LLN violation starts as early 
as A^ = 10^ — 10^ while in the valleys one has to wait until A^ = 10'^ — 10^ in order to observe 
it. (In both cases, it starts prevailing from the larger e side.) This two- fold occurrence of the 
LLN violation leads to the successive peak- valley structure around A^ ^ 10'^, which becomes 
outstanding at A^ ~ 10"^. Beyond that, up to the largest analyzed A^ (=10^), the MSD is 
independent from A^ except for e ^ 0.005 where the maps still follow LLN approximately. 

2. Mean field distributions 

Now let us discuss the mean field distributions. In the rank density plot — the bottom 
panel of Fig. — the distribution is assigned a rank as follows. 

0: The distribution is gaussianQ. The MSD is the same within 20 percent error with that 
at £ = with common a and A^. 

1: Gaussian but with a sizably enhanced MSD (the hidden coherence). The MSD can be 
even factor of ten larger than the MSD by LLN. 

2: A singly peaked distorted gaussian, or a trapezoidal distribution. 

We are interested in detecting the coherence among the evolving elements by the enhanced MSD. 
^The h{n) distribution cannot be a precise gaussian as is limited in [—xl,xl]- When we discuss 
whether it is gaussian or not, we concern whether the essence of CLT, that the convolution of 
independent random distributions peaks like gaussian, is in action or not. 
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3: Either it has a few sharp peaks on top of a broad band, or it is an apparent overlapping 
gaussian distribution. It manifestly shows the periodic motions of the maps. 

4: The distribution consists of a few sharp peaks only. 

At rank zero the h{n) distribution obeys both LLN and CLT and the maps may be thought 
as independent random numbers with a common probability distribution. Oppositely at 
rank four the maps are in periodic motion and so is the mean field. The ranks are organized 
in a way that the periodicity of the elements becomes more manifest with an increase of the 
rank. The MSD surface and the rank density plot both together reveal a simple rule: The 
MSB is high wherever the rank is high and vice versa. The rank distribution plot is almost 
a contour plot of the MSD surfac^ 

3. The regularity in the MSD enhancement 

The above rule persists for larger too. In Fig. 0, we show for N = 10^ the h{n) 
distributions at MSD peaks in the upper small boxes and at valleys in the lower. We findQ: 

(i) At any MSD peak, the rank is always high — rank three. This succinctly tells that the 
high MSD is induced by the maps evolving in quasi-periodic motion at the peak e 
values. 

(ii) On the other hand, at any MSD valleys, the rank is one (the MSD-enhanced gaussian) 
and reflects the hidden coherence. 

In short: the MSD peaks at large N come from the quasi-periodic motion of the element 
maps and the hidden coherence is restricted to the MSD valley at large N. 

^The rank assignment to each of thousands of distributions was a painful task. It was thrilling 
that independently determined two diagrams turned out in perfect match. 

^These two rules actually hold for = 10^ up to the largest A(= 10^) of our analysis. See below 
for further discussion on the N dependence of the h{n) distribution. 
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We should add that the most prominent MSD peak and the deepest valley at the front of 
it are two extremes. At the former (0.040 < e < 0.050) the distribution is either rank three 
or even four and the MSD peak starts even at small A^. The rank four distribution exhibits 
a periodic coherent motion of maps. We show below that it is due to the formation of p3c2 
cluster attractor; the lack of one cluster leads to a high MSD. At the latter, the distribution 
is also rank four but the MSD suppression is realized by the symmetrically populated p3c3 
MSCA. We will further investigate the periodicity manifestation in general shortly below 
introducing other phenomenological means too. 

4. The hidden coherence revisited 

The coherence, as is observed by the violation of LLN, occurs at any e value in the 
range 0.005 — 0.12 except that the onset of the violation is earlier at MSD peaks. [See Sec. 
III.A.l]. But the hidden coherence implies more; the MSD must be enhanced but the mean 
field distribution must remain Gaussian — the rank must be one. 

To pin down the regions of the hidden coherence on the e — N plane, let us investigate 
the change of the h{n) distribution with the increase of A^. Fig. ^ exhibits a typical case; 
the e = 0.0682 corresponding to one of MSD peaks at a = 1.90. Just when the LLN 
violation starts at A^ = 10^ ~ 10^, the rank becomes one. But notably, for A^ beyond 10^, 
the rank soon becomes two and simultaneously the MSD peak- valley structure turns out. 
For A^ ^ 10^, the rank becomes three and the peak-valley structure becomes remarkable. 
The regions of the hidden coherence are thus restricted to a very small part. Excepting the 
transitive region A^ ^ 10^ — 10'^, it has to be only MSD valleys for the h{n) distribution to 
remain Gaussian and further A^ ^ 10^ for the MSD to be enhanced. (The deepest MSD 
valley must be also excepted since we observe the apparent periodic motion of p3c3 MSCA. 

) 
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B. The periodicity manifestation in the turbulent regime 

The MSD peak-valley structure reflects a periodicity manifestation in the turbulent 
regime at various strength depending on the value of e. Let us substantiate this issue by the 
following quantities; (i) the distribution of maps and their mean fleld, (ii) map-orbits, (iii) 
the temporal correlator of mapsQ and (iv) the return map of h{n). In Fig. ^, the a is set at 
1.90 and above quantities are listed in a row for each typical iteration at characteristic e. 
The lattice size is flxed at = 10^ in order to shed light more on the predominant period 
three window than the other windows. 



1. The p3c3 MSCA: the event e = 0.036 

Let us flrst investigate the region of the deepest MSD valley. In the map-orbits we find 
clearly three lines showing the period three motion of maps in three clusters and the map 
distribution shows three delta peaks. The mean field (the black circle) is almost constant 
due to the high population symmetry and accordingly the h{n) return map shows almost 
degenerate three points. The temporal correlator oscillates in period three. All exhibit the 
formation of p3c3 MSCA. 

There is a slight subtlety that the state is actually bifurcated — six clusters of maps with 
almost equal populations in the bifurcated period three motion. This is seen by the tiny 
split in the orbits near zero|. For a = 1.90, always the bifurcated p3c3 MSCA is formed at 
e ~ 0.035, while at slightly higher e (~ 0.037 — 0.041) the final state is either a genuine p3c3 



^C{t) = (^x{n+t) ■x{n) /\x{n+t)\\x{n)\j with the relative vector x{n) = {xi{n) — hn, • ■ ■ ,xj\f{n) — 
hn)- The average (• • •) is taken over n for the last 1000 steps. 

^The six orbit points consist of three doublets of points and the two points in a doublet are very 
close each other. We have checked this numericaUy but the map distribution with the bin size 
5 X 10^^ shows only three delta peaks. 
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MSCA (90 percent) or unstable period three clusters with high rate mixing (10 percent) 
depending on the initial configurations. We come back to this point in the stability analysis 
section below. 

2. The p3c2 cluster attractor state: the event at e = 0.042 

At the nearby stronger coupling {e ^ 0.041 — 0.051), the maps almost always split into 
two clusters with the population ratio approximately 2 : 1 and the two clusters oscillate 
mutually in period three. The map-orbits sampled at e = 0.042 clearly exhibit this p3c2 
state. The mean field oscillates in period three with a large amplitude due to the lack of 
one of the MSCA and hence leads to a prominent MSD enhancement. See also the largely 
separated three points in the h{n) return plot as well as the temporal correlator in period 
three motion. Note that this high MSD is independent from the number of maps — a 
way of violating the LLN — simply because the large A^ GCML dynamics is reduced to that 
of two clusters. The MSD is solely determined by the scale given by cluster orbits and the 
population ratios. As a check let us try an estimate of the MSD. For the population ratio 
6*1 : 6*2 it is given by 

{6hy = {s/3){el + el - 1/3) - {2T/3){el + ej - 6,62), 61 + 62 = 1 (3.1) 

with S = J^^k and T = J^^k^k+i- Let us take as approximations 6i : 62 = 2 : 1 and the 
orbit points Xk at the tangent bifurcation pointQ. Then we obtain {6h)'^ = 2^/{3^7) ~ 0.169 
in good agreement with the observed value 0.16 ± 0.01. 



= 7/4 and Xk = 2/21 + 8/(3^/7)cos ((6^ + 2A;7r)/3), A: = 0,1,2 with 9 = tan-i(3\/3) for the 
stable set. Numerically (0.9983,-0.7440,0.03140). 
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3. The peripheral point to the p3c3 MSCA: two events at e = 0.032 



Here we have to account for the first transient behavior of the maps. In the event (A), 
the maps drop into p3c3 cluster attractor after a long iteration (at n ~ 8 x 10'^), while in 
(B), they remain in a few unstable clusters in mutual period three motion until the last. 
The event (A) is essentially the same with the p3c3 MSCA event. We should only note that 
the broad lower band in the h{n) distribution is an artifact of the first transient motion of 
maps. In (B), the clusters are unstable and there is a mixing of maps between the clusters; 
hence we can see only three clouds in the h{n) return map. But the mixing rate is not so 
high as we can see from a gradual exponential decayf] of the correlator with r ~ 140, which 
clearly shows a damped oscillation in period three. 

4- The variation of dynamics with e 

Let us have a bird's eye view of Fig. ^. From the row e = to e = 0.042 is the path from 
randomness to periodicity. At e = the maps evolve freely in pure randomness. We observe 
in the map distribution many sharp peaks with fractal structure. These reflect unstable fixed 
points of a single map. But the h{n) distribution — the convolution of the map distribution 
— is gaussian due to CLT. It is sharp due to LLN. The maps evolve randomly in a simple 
logistic pattern and the correlator decays almost instantly. With increasing coupling e, 
the coherence between maps is increased. The correlator reveals the precursor of the period 
three cluster attractor by its p = 3 oscillation and becomes prolonged. The map distribution 
turns into three broad bands losing sub-peaks and becomes finally sharp three delta peaks. 
Because of the increased coherence, the h{n) distribution retains the orbits structure even 
after the convolution and the rank of the h{n) distribution is gradually increased. Finally 

^We quote by r the number of steps in which the correlator decreases to 10^^ as an estimate of 
the mixing rate. 
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the rank-four distribution appears in the period three region. 

In the the period three region, we first observe the formation of the p3c3 MSCA and at 
shghtly higher e the p3c2 cluster attractor. This region continues up to e ~ 0.050. 

Beyond this, everything proceeds reversely till e ~ 0.058. The rank gradually decreases 
and the correlator gets shortened. Fig. ^ ends at this position. At larger e, the MSD shows 
small peaks and valleys in the range e ~ 0.06 — 0.1. Above s ~ 0.10 the correlator catches 
the precursor of the ordered two clustered regime. The path to the periodicity is repeated 
and eventually the period two regime starts around e ~ 0.2. 

This is the bird's eye view of the turbulent regime at a = 1.90 and N = 10^. For 
larger A^(^ 10^), the MSD surface shows peaks and valleys more remarkably. The bulk 
of above variation of dynamics with e holds also at the local scale — for each couple of 
nearby peak and valley. At the peak /i(ri)-distribution is rank three and, with the change 
of e to the nearby valley values, the rank gradually decreases down to one. At the higher 
(lower) nonlinearity a, we observe the same dynamics if the coupling is shifted to the larger 
(smaller) side with appropriate amount. For instance, the period three attractor region 
e ^ 0.032 - 0.050 at a = 1.90 is shifted to e ^ 0.06 - 0.08 at a = 2.0. This suggests curves 
of the balance in the a, e parameter space. But why the period three attractor states are 
formed at the particular e region? Aren't there any other cluster attractors with different 
periodicities? Our next task is to answer these questions deriving the curves of the balance 
analytically. 



IV. AN ANALYTIC APPROACH 

A. The tuning condition and period three clusters 

Let us consider an idealized (exact) MSCA. It is a state of GCML under three conditions, 
(i) The maps of GCML split into c clusters with an exact population symmetry, (ii) the 
synchronization of maps is perfect so that there is no variance of map positions in each of 
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the clusters, and (iii) the clusters mutually oscillate around p = c orbit points. Using this 
idealized state as a key, we derive below the tuning condition for the MSCA formation. For 
brevity we explain our approach with respect to the p3c3 MSCA in detail but everything 
below also goes through for the other MSCA with higher periodicity. In a p3c3 MSCA three 
clusters A, B, C move cyclically round three fixed positions Xi, X2, X^. Such a system of 
orbit points exists as a triple intersection point of three surfaces given by 



S,; : Xi 



, (2,j,A;)G {(1,2,3), (2,3,1), (3,1,2)}. (4.1) 



At a = 1.90, e = 0.040, for instance, we have two solutions (0.96301,-0.00499,-0.72851) 
and (0.95521, 0.07076, —0.69993); the former is stable and the latter is unstable. In such an 
exact MSCA, the meanfield h{n) is a time-independent constant: 

1 ^ 1 1 ^ 1 3 

Mn) = ^E/(^^H) = o E /(x,(n)) = -E/TO = TTE^^ = ^% (4-2) 

1=1 I=A,B,C i=l i=l 

where Xj{n) denotes the coordinate of the cluster / at time n and the last equality follows 
from ( ^.21 ) or ( |4.1| ). Therefore, if MSCA is produced, the GCML evolution equation ( ^.Ij ) 
becomes 

x,{n + l) = {l-E)faix,in))+eh* (2 = 1, ■ ■ ■ , AT), (4.3) 
fa{x) = 1 - ax'^, 

where the time-dependent term h{n) is replaced by a constant h*. Every one of the maps 
evolves by a common equation at each step in ( p.l|) and further by a unique constant equation 
in ([4.3|) . As is noted by Perez and Cerdeira some years ago, we can cast this unique 



equation to a standard logistic map with a reduced nonlinear parameter b 

y,{n + l) = l-b {y,{n)f (z = 1, ■ ■ ■ , AT). (4.4) 
by a linear scale transformation 

Viin) = {l-6 + eh*)-^Xi{n), (4.5) 
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and the reduction rate of the nonhnearity parameter is given by 

r = - = (l-e)(l-e(l-h*)). (4.6) 
a 

At MSCA, the mean field h* is constant, so the reduction factor r is also constant. If 
the clusters of MSCA oscillate in period three, so do the maps yi{n) — the two solutions 
(xi, X2, Xs)^'^\ {u = 1, 2) of the cyclic equation ( [4.1|) agree with the two sets of period three 
orbit points {y{n),y{n + l),y(n + 2)Y'^\ (z/ = 1,2) of the logistic map (|4^ ) modulo the scale 
factor in ( [4.5|) . The reduction factor r must reduce the high nonhnearity a of GCML down 
to the b in the period three window. It starts at 6th = 7/4 by the tangent bifurcation and, 
after sequential bifurcations and windows in the window, it closes at 6 = 1.79035 by the 
crisis. The range of the period three window 6 = 1.75 — 1.79035 requires a reduction factor 
r in the range 0.942 — 0.921 for a = 1.90. Each r within this range gives a constraint curve^ 
on e,h* plane via ( |4.6| ). 

There is another constraint from self-consistency; the average value y* of the transformed 
maps must also obey ( [4.5| ) so that 

y*^}-j2y, = {l-e + eh*)-'h*. (4.7) 
i=i 

Here y* is a function of the nonlinear parameter b — it is simply an equal weight average of 
the period three stable orbits of the single logistic map ( [4.4|) and can be estimated solely by 
the property of the logistic map without any recourse to the GCML evolution equation. At 
a given y* this again gives a constraint curve on e, h* plane. Let us work out the e at the 
intersection of the two curves. By eliminating h* from ( [4.6| ) and ( [4.7| ) we obtain 

^ = 1 - ^ - \/r(i-r) + (^)'. (4.8) 

and both r = b/a and y* in the right hand side are determined by b. This is the tuning 
condition. This predicts the necessary value of the coupling e for the GCML at a given a to 



^It is possible to transform formally ( |2.1| ) to a standard form at each step but then the reduction 
factor r may fluctuate step by step. Then it does not single out a line. 

16 



form MSCA due to the periodic attractor of the single logistic map with b. The function y*{b) 
is a well-known square-well, y* ~ 0.284 at 6 = 1.735 slightly below the tangent bifurcation 
point 6th, and it drops sharply {y* — y*^ oc y/bth — h) at 6th- From matching of the coefficients 
in {hf{y) -y = b\h{y) - y) {liUiv - yi))\ we obtain = 1/(3 ■ 26) = 2/21 = 0.095. 
Similarly y*{b) = (1 — v46 — 7)/66 up to the first bifurcation point 6 = 1.769. Then, y*{b) 
varies smoothly|3 around 0.08 until the end of the window (6 = 1.7903) and finally increases 
sharply {y* ^ 0.18 at 6 ~ 1.793). This y*{b) put into (|4.8|) gives the following estimates of 
£ for a = 1.90; 



A : 


e 


= 0.0514 


at 




= (1.735, 


0.284), 


r 


= 0.913 




B : 


e 


= 0.0422 


at 




= (1.750, 


0.095), 


r 


= 0.921 




C : 


e 


= 0.0363 


at 


{b,y*) 


= (1.769, 


0.069), 


r 


= 0.931 


(4.9) 


D : 


e 


= 0.0305 


at 


{b,y*) 


= (1.790, 


0.080), 


r 


= 0.942 





The estimates A, B, C and D are respectively below the threshold, at the threshold, at the 
first bifurcation point in the window, and at the closing point of the window. Note that the 
route A ^ D is in the direction of increasing 6, which in turn is the direction of decreasing 
coupling constant e, since the larger 6 requires only a smaller nonlinearity reduction. We 
should stress that the tuning condition ( |4.8| ) is a necessary condition. For the p3c3 MSCA to 
be stable, the orbit of the reduced logistic map must be also stable]^. Period three logistic 
orbit still continues to exist even beyond the first bifurcation point C, but it is unstable. 
Therefore, for an exact p3c3 MSCA to be formed, the e range must be within the estimates 

^^The largest rapid variation is the tiny anti square well (Ay* ~ 0.01) due to the 3x3 window 
at 6 = 1.7858 - 1.7865 

^^As we show in the stability section below, the Lyapunov exponents of GCML at the p3c3 MSCA 
consist of A — 3-fold degenerate one and three in general non-degenerate ones. For the p3c3 MSCA 
to be stable, at least the former degenerate exponent must be negative, which implies the reduced 
map orbit must be stable. That all these exponents are negative at MSCA is shown also below. 
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C-B, but neither within D-C nor beyond D. Similarly an exact bifurcated p3c3 MSCA must 
be formed within D-C. Our tuning condition does not guarantee the formation of the MSCA 
but it does limit the £— range in which the formation is possible. The observed ranges of the 
p3 cluster attractors are listed in Table |. At a = 1.90, p3c3 MSCA is formed in the range 
e ~ 0.037 — 0.041 and its bifurcated state in e ^ 0.032 — 0.037. The predicted ranges are 
respectively e = 0.0363 — 0.0422 and e = 0.0305 — 0.0363. In both cases, the agreement is 
remarkable and we see that the formation actually occurs at any allowed e value. 

As for the p3c2 cluster state, we need a caution in using the tuning condition. It is derived 
under the assumption of the constancy of the mean field. Thus, as a matter of principle, 
it cannot be applied for the asymmetrically populated state. However, the p3c2 state is 
formed with a slightly higher coupling e and the orbits of two clusters are approximately the 
same with the MSCA orbits. Therefore, the p3c2 cluster attractor is certainly still under 
the control of the period three window. We estimate the range by the extension of the 
period three window at the higher coupling side B-A — the intermittency region. This gives 
e = 0.0422 — 0.0514 in good agreement with the observed range of the p3c2 cluster attractor 

0.041 — 0.050). It is interesting to note that the GCML final states at this e range 
actually consist of two types depending on the initial condition; the p3c2 cluster attractor 
(~ 80%) as well as the unstable period three clusters with mixing of maps (the rest). See 
Fig. P below. The estimate by B-A relates intriguingly the intermittency of the element map 
to the GCML phase of co-existent stable and unstable periodic clusters. We are aware that 
we cannot take the success of the estimate for p > c states on the same footing with that 
for the MSCA but at least it gives a good rule of thumb for the p > c state. 

B. Foliation of the logistic windows in the turbulent regime 

What is the case for other a values? Do the other windows also show up in the expected 
e range in the turbulent regime? To check these systematically, let us note that the tuning 
condition defines a one-parameter (b) family of curves in the model parameter space (the 
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a, e— plane) of the GCML. Each curve is labeled by b and written as a function of the 
reduction factor r as 



(a('')(r),£(')(r)) 



b ^ _ ry*{b) 
r' 2 



V 



r{l-y*{b))+ ('^^^] I, r<l. (4.10) 



It emanates from the point {a^^\ e^^'^)\r=i = 0) and with the decrease of r it develops in 
the parameter space in the direction in which both a and e increase in a certain balance. 
If our success above is a general one, all of the GCML with the parameters being set at 
(r),e(^)(r)) along a curve labeled by b should be commonly controlled by the same 
dynamics of the single logistic map at b. 

In Fig. 1^ we find that this is indeed the case. Each panel shows the MSD of the h{n) 
distribution as a function of e at a given a as well as the expected zones for the manifestation 
of the outstanding six windows in Table.|I|. The curves ( [4.10|) are displayed underneath the 
panels and link the respective zones. At each zone, a MSD valley due to MSCA should 
appear in the lower e side and a MSD peak hj p > c cluster attractors at the nearby higher 
e. We find that this works with almost no failure in all panels and with respect to all six 
windows. 

The effects of the logistic windows propagate along the curves ( ^.1CI| ), which may be 



called as foliation curves. The curve with the label b links together those GCML commonly 
subject to the same logistic window dynamics at b. Accordingly the family of the curves 
produces the foliation of the single map dynamics. The foliation occurs because, under the 
global interaction, the maps of the GCML form a macroscopically coherent state. Even 
though the coupling in the turbulent regime is very small, the coherence prevails over the 
GCML maps if the tuning condition is met. A few remarks are in order. 
1. Desynchronization along the foliation curve The periodicity manifestation becomes 
weakened at a higher reduction and there is a threshold rth ~ 0.95. For r ^ 1, both MSCA 
and associated p > c clusters are formed in tight synchronization. Towards rth ~ 0.95 the 
clusters broaden. There is no mixing yet among the clusters but the maps move chaotically 
in each cluster. Below rth, we observe only the periodicity remnants — on one hand the 
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overlapping-gaussian /i(n)-distribution along the curve which had the p > c cluster above 
rth, and on the other hand the MSD-enhanced gaussian (the hidden coherence) along that 
of MSCA. See Fig. |. 

2. Left-right asymmetry of the MSD curves The MSD curves in Fig. ^ (and |^) show an 
interesting feature — in each panel the smaller e region (the left ) has an ample amount of 
peaks and valleys, while the larger only a few broad ones. As for the single logistic map, on 
the other hand, there are as many windows in the smaller b as in the larger. 

This is naturally understood by the difference in the reduction factor r between the zones 
in a panel. In a way, each panel is a screen which displays the windows of the single logistic 
map by using a macroscopic coherent state of GCML. But the panels set at fixed a values are 
inclined — a smaller e (the left) implies less reduction, i.e. r ^ 1. The rth divides the panel 
at e ~ 0.030 via ([4.8|) . The left sensitively displays the sharp peak- valley structure induced 
by cluster attractors. The right, on the other hand, can reflect only the accumulation of 
the periodicity remnants from nearby windows, being dominated by the prominent one at 
its respective zone. As a check we set the panels at fixed r values. Then they displayed 
windows without asymmetry, and with a higher sensitivity at r closer to one [?]. 

3. Cluster attractors with higher periodicity Let us search cluster attractors with periodicity 



higher than three. Here we give two samples in Table.[lTl. 

p = 4 clusters These appear in the left most zone in the a = 1.95 panel. From the 
window b data in Table.|T| the necessary reduction from a is very small — r ^ 0.995 — 
so we expect definite clusters. We indeed find the expected sequence of clusters^ p = 4, 
c = 4(MSCA)^3^ 2 in tight synchronization at the right e. 

p = 5 clusters There are two p = 5 windows in Table.||. We choose the one at the lower 
b and set b = 1.66 which amounts to r = 0.980. Since r is in the mid of one and rth, we 



^^The single cluster cannot be formed. The focusing by averaging does not act there and the tiny 
variance is instantly amplified. It appears far in the coherent phase (e ^ 0.4 at a = 1.90). 



20 



expect the clusters are not in complete synchronization but yet there is no mixing of maps. 
Indeed the sequence of attractors p = 5, c = 5(MSCA) ^ 4 — > 3 is observed at the expected 
e and it terminates before the lowest one (j95c2). 



V. STABILITY OF THE PERIOD THREE CLUSTERED MAP STATES 

Here we adopt the Lyapunov analysis. As one superlative ability, it can be applied to 
both diverging and converging system orbits so that it can detect the possible coexistence 
of multifold finial states depending on the initial configurations. We measure the maximum 



Lyapunov exponent Amax by a standard method |18| which keeps track of an A^-dimensional 



shift vector (5x(n) evolving under the non-autonomous linearized equation associated with 



5xi{n + 1) = -2a I ^ ^ + ^) Xi{n)Sxi{n) + ^ X! ^iH'^^i(^) | • (^-l) 

The Amax is the average of the logarithm of the expansion rate of the shift vector (with 
intermediate renormalizations). For both Amax and MSD, we discard the first transient 10^ 
steps — generally it takes only 10^ ~ 10^ steps for the cluster formation. 

Let us first check the e— dependence of the stability of attractors. We choose = 10^, 
fix a at L90, and vary e in the range 0.030 — 0.052 with the inclement Ae = 10~^. We show 
in Fig. I the Amax in the upper and the MSD in the lower^. 

We observe in the Amax plot three remarkable structures of low Amax events, 
(i) A seagull structure {e = 0.032 — 0.037) with a sharp cusp at e = 0.0352 — all events are 
bifurcated MSCA with good population symmetry {Nj/N (1 ± 0.05)/6). Note that the 
events form also a seagul in the MSD and the cusp positions agree precisely. The bifurcated 



For reference, the CPU time for N = 10 GCML is approximately two minutes for one mea- 
surement of Amax (2^^ steps for precision 10^^) plus MSD (10^ steps) on a modest supercomputer 
VPP300/6. The total is 2 min x 40(initial configurations) x 220 (e-values) 300 hours. 
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MSCA is the more stable if the mean field fluctuation is the less and it is the most linearly 
stable (Xmax = —0.38^ with the minimum fluctuation (5h^ 2 x 10~^/ 

(ii) The first low band {e = 0.037 — 0.041) — The p3c3 states. The population distributes 
around the exact MSCA — Oi = Nj/N (1 ±0.15)/3. The events near the lower boundary 
(Amax < 0) are p3c3 events with good population symmetry and with low MSD. 

(iii) The second low band {e = 0.041—0.051) — the j93c2 cluster attractor. The corresponding 
MSD is, contrary to (ii), extremely high because of a lack of one cluster to minimize h{n) 
fluctuation [See Sec. III.B.2]. 

The foliation of the critical points A, ■ ■ -,0 from the period three window defines three 
e-regions I(D-C), II(C-B), III(B-A) [Eq. ( |4.9| )]. The region I is the allowed region for the 
formation of the bifurcated MSCA (MSCA*), II the p3c3 MSCA and the p3c2 attractor 
cluster is expected in III. As we see clearly in Fig. H, the regions I, II and III respectively 
embody the structures (i), (ii) and (iii) just in agreement with our prediction. 

Let us note a remarkable feature in the events in the two wings of the seagull (i). Here all 
events come out with positive Amax (~ 0.1 — 0.2). For a system with low degrees of freedom, 
the positive Amax implies chaos. But here even with positive Amax, the maps always form 
bifurcated p3c3 state. There is actually no contradiction. The global motion of the clusters 
is periodic, but, inside each cluster, maps are here evolving randomly with tiny amplitudes 
(^ 10~^) in sharp contrast against the complete synchronization at the cusp. The Lyapunov 
exponent measures the linear stability of the system with respect to the small deviation of 
the element position. It is sensitive to the microscopic motion of the element of the system 
and hence yields the positive exponent. But for a larger deviation, nonlinear terms can 
become relevant and pull back the mapQ. This type of map motion — microscopically 
chaotic but macroscopically in the periodic clusters — may be called as confined chaos. 



We have verified this by inputting pulses on randomly selected maps. The analytic formulation 
of the nonlinear effect is most wanted for. 
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We hereafter devote ourselves into the investigation of two outstanding structures, 
namely the bifurcated MSCA seagull and the p3c2 cluster attractor. 

The linear stability analysis of the bifurcated MSCA In order to understand the salient 
cusp at £ = 0.0352 in the Lyapunov exponent plot, let us consider the linear stability matrix 
of the GCML. 

(1) For the configuration of maps in six clusters, the N x N linear stability matrix of the 
GCML for evolution of one step can be written as 



Ml = (!-£) 



XiEi 









■ • • XqEq 



XiHii ■ ■ ■ XqHiq ^ 



+ 



N 



X^H, 
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XrH; 



6-" 26 



XiHqi ■ ■ ■ XqHqq j 



(5.2) 



multiplied by an overall factor —2a, where Xi{I — l,---,6) are the coordinates of the 
clusters, Ej aNiX Nj unit matrix, and Hjj is a A^/ x Nj matrix with all elements one. The N 
eigenvalues of Mi consist of two sets. One is a set of 6 eigenvalues A*^'^'' = —2a{l — e)Xi, (/ = 
l,---,6), each (Nj — l)-fold degenerate. The degenerate eigenvectors of A^^^ are of the 
form Col.(0; • • • ; 0; (1, 0, • • • , 0, —1, 0, • • • , 0); 0; • • • ; 0), that is, all column blocks, each for 
one cluster, are fulfilled by except for the I-th block which has 1 as the first element 
and —1 as one of the other Nj — 1 elements. The eigenvector of A^^^ represents a shift of 
the system orbits within the I-th cluster. The other is a set of (in general non-degenerate) 
6 eigenvalues A/, which are the same with the ones of the 6x6 stability matrix Mi; red 
associated with the cluster evolution 



Xj{n +!) = (!- e)f{Xj) + s ^ OjfiXj), (/ = 1, • • • , 6). 

j=i 



(5.3) 



The Ml; red ioT the cluster dynamics is derived from Mi hj Ej 1 and Hjj — > 6j. The 
eigenvector of Mi subject to A/ is (^f 1; • • • ; Ce^)) "^i^^ (Ci; ' ' ' > Ce) being that of Mi; red- 
(2) The stability matrix Mp of GCML for the evolution of p steps is given by the chain 
product of p of Ml along the system orbit. The eigenvalues of Mp again consist of two sets. 
One is the set of 6 eigenvalues A^^^ — [— 2a(l — £)]^ 11^=1 Xj, each (Nj — l)-fold degenerate. 
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and the first set eigenvectors of Mi remain tlie eigenvectors of this set. The other is the 
same with the ones of the Mp.j-ed — the p-th iterate of Mi-red- This mechanism holds at any 
population composition among the GCML clusters. 

(3) Now, when the population symmetry among the clusters are exact, all of the maps obey 
a unique quadratic mapping ( |4.3| ) with a constant mean field h* which is equivalent to a 
standard logistic map ([4.4|) with a reduced nonlinearity b via the scale transformation ( [4.51) . 
For b from the first to the second bifurcation point in the p = 3 window(66 = 1.76852915 
to bi2 = 1.777221618), the reduced map y evolves in period six and so do the GCML six 
clusters. This is the bifurcated MSCA. We can write the correspondence as 



yo 



yi = fbivo) 



y, = {hf{yo) ^ yo = iUYiyo) 



X, 



\X,j 



Xi 



\X,j 



Xo 



\Xej 



(5.4) 



In the MSCA, all of the six eigenvalues of Mg in the first set degenerate into a single value 
A = [-2a(l - €)fUhiXi with degeneracy ELi(^/ - 1) = N - 6. By (0), (p|) and 
( [4.6|) we find A = (—26)^ 11^=1 Z/i; that is, A is nothing but the Lyapunov eigenvalue of the 
single logistic map for the p = 6 motion. As for the other set, the Mq. for the symmetric 
configuration 9j = 1/6 is a chain product of six matrices, that is, Mf. red^i- red ' ' ' ^i- red 
with 



Ml 



red 



-2a 



V 



'l~e + r])Xi 7]X2 ••• r/Xg 
r/Xi (l_£ + ry)X2 ■■■ r/Xg 

vXi {l-s + 'n)X^ 



\ 



^=6 



(5.5) 



and other five matrices are obtained by cyclically changing the orbit points Xj by ( |5.4|) . By 
a simple algebra using ( [4.5|) and (|5.4|) , we find that the eigenvalues of Mg in the second set, 
which are in turn the ones of Mg; red, are A with corrections of order rj. 
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(4) Now we are ready to work out the seagul cusp position. Because the A is proportional to 
the product of the period six orbit points of the single logistic map fb, it becomes zero when 
one of the orbit points becomes zero. At this very instance, the — 6 Lyapunov exponents 
become —oo and the other 6 exponents become also very small proportionally to log(|r7|)/6. 
The 6 is a solution of (0) = and the relevant solution be = 1.772892 gives Ec = 0.035192 
and Amax = —0.361519 for a = 1.90 — both are in remarkable agreement with the observed 
cusp of the GCML Lyapunov exponent. 

Over the seagul e— range, Mg; red has four complex [{\k,\l),k = 1,2)] and two real 
eigenvalues and gives four exponents. The Amax is given by one of the two sets of complex 
eigenvalues, while the (A^ — 6)-fold degenerate exponent from A runs in the mid of the four. 
The predicted Amax is shown in Fig. ^ and explains the data well. The slight deviation off 
the cusp is due the small population unbalance; it is the larger for the larger MSD events. 

The dependence of the Xmax on the population ratios We proceed with the following 
algorithm after detecting the clusters by the gaps. The six MSCA* clusters evolve in the 
bifurcated orbits of p3c3 MSCA. They can be regarded as three doublets — (C/i, Ch), ^ = 
1,2,3 so that the two clusters Cj^ and Cj^ in a doublet evolve close together. We combine 
the two populations in a doublet into one and define s, t and u as {Nf^ + Nj^)/N in the 
decreasing order. 

In Fig. ^(a), we exhibit the averaged —Amax on the s, t— plane from the 2 x 10^ random 
events for A^ = 10^ GCML with a = 1.90, e = 0.035. At the top of the pyramid-shaped 
surface the Amax is negative and at its minimum. It occurs precisely at the most symmetric 
population configuration and we find only an event with almost perfect population symmetry 
is formed. The Amax is negative over the bulk of events around the symmetric point — MSCA 
is linearly stable. The exception occurs only near the boundary (the round curve), where 
the Amax is mostly positive and small (Amax ~ 0.05) and the maps form the confined chaos. 

The p3c2 Cluster Attractor We have done a similar high statistics analysis at a = 1.90, 
e = 0.048 for the same A^ = 10^ GCML. The final states are two fold; p3c2 cluster attractor 
(83%) and the unstable p3 clusters with mixing (the rest). Hereafter we analyze the former 
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in Fig. ^. In the region 0.55 < 9 < 0.61, the p3c2 clusters are tightly bounded and linearly 
stable. Here the dynamics of the GCML is reduced to that of two clusters. Just like the 
p2c2 state in the ordered two clustered phase , the p3c2 orbits bifurcate with the change 
of — the ratio 9 can be used as a control parameter even in the turbulent regime. 

However, there is a remarkable difference too. In p2c2 there is no stable attractor for 
9 outside the window. In the turbulent regime, on the other hand, a loosely bound p'ic2 
state can be formed — the three orbit bands in the edge regions. This state is again the 
confined chaos. The Amax is positive (0 < Amax < 0.2) and the maps fluctuate randomly 
in each of the two clusters. But the clusters are in a macroscopic period three motion. As 
the probability distribution shows, this is formed as frequently at the p3c2 cluster attractor. 
The state of confined chaos at the unbalanced population is a characteristic feature of the 
cluster attractors in the turbulent regime. 

VI. CONCLUSION 

In this article we have revisited the GCML of the logistic maps and studied in detail 
its so-called turbulent regime. We have presented our new phenomenological findings in an 
extensive statistical analysis, which as a whole tell that the turbulent regime is under the 
systematic control of the periodic windows of the element logistic map. In particular we 
have shown that the hidden coherence occurs only in a very limited regions in the turbulent 
regime. 

There appears remarkable p3c3 MSCA states as well as p3c2 cluster attractors induced 
by the period three window of the element map. Our tuning condition predicts by a family of 
curves how the dynamics of the element map foliates in the parameter space of the GCML. It 
successfully explains the salient peak-valley structures of the MSD surface and tells us where 
to see the remarkable sequence of the cluster attractors of the type p,c = p —>■ {p — 1) — >■ 
(p-2)^---. 

We have also investigated the linear stabihty of the period three cluster attractors. Both 
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the p3c3 MSCA and its bifurcated state are linearly stable when the population symmetry 
is good and MSD of the meanfield is minimized. We have analytically explained the value 
of the coupling £ at a given a for the formation of the most stable bifurcated MSCA. The 
p3c2 cluster attractor is also linearly stable in the 6'— window even though the MSD of 
the h{n) is quite high. For the unbalanced population configuration the system forms an 
interesting state of confined chaos which is a characteristic feature of the cluster attractors 
in the turbulent regime. 

There remain interesting unsolved problems. One concerns with the state of confined 
chaos newly found in the turbulent regime. It is a state consisting of a few clusters in 
macroscopic periodic motion and maps move around chaotically inside each clusters. Re- 
garding the linear stability the Lyapunov exponent is positive. It is tempting to single out 
the nonlinear effect which confines the maps in periodic clusters. A related problem is the 
onset of the incomplete synchronization with the decrease of the reduction factor r along the 
foliation curves. The other concerns with the variation of the dynamics with the system size 
N. We have found that the system becomes an extremely sensitive mirror of the element 
dynamics with increasing N. The salient evidences are shown in Fig. ^ and Fig. ^ but we 
are unable to explain why so. In field theory the vacuum at the spontaneous break down of 
the symmetry is stable only when the degree of dynamical degree of freedom is infinite []19| . 
If we may regard the randomness of GCML maps as a symmetry, the MSCA with no h{n) 
fluctuation corresponds to a vacuum at the symmetry breakdown and the formation of it by 
synchronization the onset of the ordered parameter. The resolution of the finite size effect 
in GCML is so tempting since it may bridge the synchronization of the maps and onset of 
the order parameter in the field theory in quantitative terms. 

As a whole this work is an exploration of order in the chaos and we have found that the 
turbulent regime of GCML is controlled by the foliation of the single logistic dynamics. 
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FIGURES 

FIG. 1. The MSD surface of mean field fluctuation (top) and the gray-scale density plot for 
the rank of distributions (bottom) on the e — N grid, a = 1.90 and the inclement As = 2 x 10~^. 
The rank varies from zero (black) to rank four (white). 

FIG. 2. Central diagram: MSD at a = 1.90 as a function of £ over e = — 0.10 with inclement 
10-^. (a) N = 10^, (b) 10=^, (c) 10^, (d)10^ (e) 10^. Top boxes: The h{n) distributions at MSD 
peaks. All are rank three. Bottom: the same at valleys. All are rank one and corresponds to the 
hidden coherence. 

FIG. 3. The variation of the h{n) distribution with the system size N ai a = 1.90, e = 0.0682 
and N = 3 X 10^,5 x 10^,2 x 10^, 10^ from left to right. Each is sampled in 10^ iterations from 
random start discarding the first 10^ steps and compared with a reference distribution (rank 0). 
The rank is given on top. 

FIG. 4. The variation of GCML dynamics with e through the prominent MSD peak region. 
N = 10^ GCML with a = 1.90 and total iteration steps are 10^ for each run. (a) The mean field 
distribution (marked as h) sampled discarding the flrst 10^ transient steps and the map distribution 
(x) averaged over the last 2 x 10'^ steps, (b) Clustering pattern. The lines are orbits of randomly 
selected 10^ maps for the last seven steps and the black circle is the mean field, (c) Temporal 
correlator between the two relative-coordinate vectors of maps. Averaged over the last 10^ steps, 
r « 0, 30, 70, 140, oo, oo, oo, 50, for e = 0, • • • , 0.056. (d) The return maps for the same steps as 
(c). The arrows indicate the period three clusters. 

FIG. 5. The MSD of the mean field distribution plotted as a function £ for = 10^ GCML 
in panels at a = 1.8(back), 1.85, ••• ,2 (front), e = — 0.12 with inclement 10~^. The foliation 

curves predicted for p = 4,5, 3, 7, 5 windows flow underneath panels and link the respective shaded 
foliation zones. For each zone a MSD-valley at lower e and a peak at higher. 
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FIG. 6. The Lyapunov exponents (upper) and the MSD (lower) measured for 40 random initial 
configurations at each e between 0.0300 and 0.0520 with inclement 10~^. a = 1.90 and N = 10^. 
The events in the seagull structures are in the bifurcated MSCA states, the first lower band events 
the MSCA, and the second ones the p3c2 attractor. These are formed in regions I (D-C), II (C-B) 
and III (B-A) separated by dashed lines just as expected. In bands M, the p = 3 clusters are 
unstable and maps are mixing among them with r « 100. The dashed line; the predicted Amax for 
the exactly symmetric MSCA. Arrows; the predicted positions of the most linearly stable MSCA 
states. 

FIG. 7. (a) The maximum Lyapunov exponent of the bifurcated MSCA over the s — t plane. 
Amax for events in one st— bin are averaged, the sign is reversed (-Amax) and the part Amax > is 
truncated for a bird 's eye view, a = 1.90, e = 0.035, = 10"^ and totally 2 x 10^ events, (b) The 
triangle APMQ is the allowed (s,t) region by the constraint s>t, t> {1 — s)/2, 1 > s + t. The 
bifurcated MSCA events accumulate in the top of the shaded tiny AAMB, which is again a tiny 
part of APMQ. 

FIG. 8. (a) The orbits of p3c2 cluster attractor for the last 512 steps, (b) the maximum 
Lyapunov exponents, (c) the number of events in each bin {A9 = 10~^) plotted at the population 
ratio 9 for 18341 events (83%) with the stable cluster formation among totally 22000 events. 
N = 10'^, a = 1.90, e = 0.048. Dotted lines to draw the eye and separate the tightly bound cluster 
states (central) and confined chaos (two edges). 
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TABLES 

TABLE I. The p = 3 periodicity manifestation in the turbulent regime, a = 1.90, N = 10^-10^. 



e— range 



Prediction 



MSDsurface 



state 



0.032 - 0.0352* - 0.037 
0.037 - 0.041 
0.037 - 0.041 
0.041 - 0.050 



0.0305-0.035192-0.0363 
0.0363-0.0422 

0.0422-0.0514 



deep vahey 
lower band 
upper band 
prominent peak 



bifurcated p3c3 MSCA 

p3c3 MSCA 
p3c2 cluster attractor 
p3c2 cluster attractor 



*The downwards cusp position of the MSD valley. 



TABLE II. The outstanding windows and their b, y*{b) values. 



period 


A: inter mittency * 


B: lower threshold 


C: the first bifurcation 


D: closing point 


Width'^ 


7 


1.5740, 0.3943 


1.5748, 0.3857 


1.5754, 0.3846 


1.5762, 0.3847 


0.0014 


5 


1.6220, 0.3610 


1.6244, 0.3077 


1.6284, 0.3012 


1.6333, 0.3032 


0.0089 


7 


1.6735, 0.3189 


1.6740, 0.2676 


1.6744, 0.2678 


1.6749, 0.2677 


0.0009 


3 


1.7350, 0.2836 


1.7500, 0.0952 


1.7685, 0.0685 


1.7903, 0.0800 


0.0403 


5 


1.8597, 0.1823 


1.8606, 0.0984 


1.8614, 0.0990 


1.8623, 0.0987 


0.0017 


4 


1.9390, 0.1287 


1.9406, -0.1633 


1.9415, -0.1668 


1.9427, -0.1657 


0.0021 



*Thc starting point of intermittency. 
"^Ab = bB-bD 
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TABLE III. Samples of the sequences of attractors with p > c. N = 10^ and a in parentheses. 
p = A: b: 1.9406 - 1.9427^ ( a = 1.95, r ^ 0.995) 

3 2 
0.0022 - 0.0024 0.0024 - 0.0026 

0.0022 0.0030 

p = 5: b: 1.6244 - 1.6333^ ( a = 1.66, r 0.980 ) 

4 3 
0.0114 - 0.0124 0.0124 -0.0130 

0.0112 0.0140 

'^The range of the window from B to D, i.e. intermittent region not included. 

'^The e— range predicted by the tuning condition and the window data in Table |I|. The 

predictions are from C-B for MSCA and B-A for p > c states. 

■^For most events, the attractor is either bifurcated or consists of five bands. 

'^D-C for MSCA and C-A for p > c cluster attractor taking into account the bifurcation. 



c 4 (MSCA) 

s 0.0019 - 0.0022 

Epred" 0.0020 - 0.0022 



c 5 (MSCA)= 

s 0.00986 - 0.0118 

ffp^e/ 0.00950 - 0.0112 
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